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The clock-like properties of pulsars moving in the gravitational fields of their 
unseen neutron-star companions have allowed unique tests of general rela- 
tivity and provided evidence for gravitational radiation. We report here the 
detection of the 2.8-sec pulsar J0737— 3039B as the companion to the 23-ms 
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pulsar J0737— 3039A in a highly-relativistic double-neutron-star system, al- 
lowing unprecedented tests of fundamental gravitational physics. We observe 
a short eclipse of J0737-3039A by J0737-3039B and orbital modulation of 
the flux density and pulse shape of J0737— 3039B, probably due to the influ- 
ence of J0737— 3039A's energy flux upon its magnetosphere. These effects will 
allow us to probe magneto-ionic properties of a pulsar magnetosphere. 

Double-neutron-star (DNS) binaries are rare, and only six such systems are known. How- 
ever, they can provide wonderful laboratories for the study of relativistic gravity and gravita- 
tional radiation (0). The discovery of such systems has been a prime objective of pulsar surveys 
since the first, PSR B 1913+16, was discovered 29 years ago ©. Apart from their rarity, they are 
particularly difficult to find because of the large orbital acceleration experienced by the mem- 
bers of the systems, resulting in large and varying Doppler effects in the observed rotational 
period. 

The pulsar J0737— 3039 was recently discovered as part of a high-latitude multibeam survey 
of the southern sky using the 64-m Parkes radio telescope ©. It was found to be in a 2.4-hr 
eccentric orbit with another compact object that the observed orbital parameters suggested was 
another neutron star. The short orbital period and compactness of the system and the high 
timing precision made possible by the large flux density and narrow pulse features of this pulsar 
promise to make this system a superb laboratory for the investigation of relativistic astrophysics. 
Relativistic advance of the angle of periastron has already been measured and orbital decay due 
to gravitational wave emission will be measured with high precision within a few months. The 
resulting in-spiral will end in coalescence of the two stars in about 85 My ©. This discovery 
significantly increases the estimates of the detection rate of DNS in-spirals by gravitational 
wave detectors (0121). 
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The discovery of PSR J0737— 3039B. No other pulsar was detected in the 4.5-min discovery 
observation of PSR J0737— 3039. However, analysis of data acquired subsequently for detailed 
studies of the system has revealed the occasional presence of pulsations with a period of 2.8 
seconds. This pulsar, henceforth called PSR J0737— 3039B (or simply B in this paper), has the 
same dispersion measure (the integrated column density of free electrons along the line of sight) 
as the original pulsar, henceforth called PSR J0737— 3039A (or simply A) and shows Doppler 
variations in the pulse period which identify it as the companion. 

The data obtained by Burgay et al. © for their studies of A used the Parkes radio tele- 
scope at a frequency of 1390 MHz and have now all been reprocessed to enable investiga- 
tion of the properties of B. Additionally, the system was observed on seven days in November 
2003 simultaneously with dual-polarization receivers centered on frequencies of 680 MHz and 
3030 MHz ©. Observations were also made of both pulsars at 1396 MHz with the 76-m Lovell 
telescope at Jodrell Bank, also with a dual-polarization receiver. Timing observations of both 
pulsars have been made covering the 7 months from May to November 2003, with individual 
data spans between 10 minutes and 5 hours. 

The main features of the observations of B are shown in Fig. 1, which illustrates the strength 
of the 2.8-sec pulsed emission in the three frequency bands as a function of orbital phase and 
pulse phase over the full orbit. Most notable is the variation in the received flux density of the 
pulsar, that is clearly visible for two brief periods of about 10 minutes duration each centered 
upon orbital longitudes 210° and 280°. Within these bursts, the strength is such that most 
pulses are detected individually. In addition to the orbital phases where it is always detected, 
the pulsar often shows weak emission elsewhere, notably between orbital phases 0° and 20° 
(Fig. IB). It is likely that observations with greater sensitivity will reveal more widespread 
emission throughout the orbit. The pattern of the visibility is essentially stable from orbit to orbit 
and from band to band over the full frequency range of 680 MHz to 3030 MHz. The absence 
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of pulsar B in the discovery observation of A can now be understood, as that observation was 
made at a longitude of 146°, where there is little emission from B. 

There are significant changes in the shape of the pulses of B with orbital phase (Fig. 1). The 
data at 1390 MHz show that it varies from a narrow intense main pulse with a weak precursor in 
the burst around longitude 210°, to a roughly-equal double-component profile with a somewhat 
greater separation near longitude 280°. It becomes essentially a single pulse at around longitude 
0°. There is also an indication for frequency evolution of pulse shape within the bursts. The 
double nature of the profiles becomes more evident at the higher frequency of 3030 MHz and 
less so at 680 MHz, where the components seem to have nearly merged. 

In order to make timing measurements of the pulsar, a phase-dependent set of pulse tem- 
plates was created by forming average profiles for each of the three burst regions. These tem- 
plates were aligned in such a way that the trailing, usually more intense, components were at the 
same pulse phase. Cross -correlation of data with appropriate templates from this set was car- 
ried out to produce times of arrival (TOAs). These TOAs were then compared with a model for 
the astrometric, rotational and orbital parameters using the TEMPO pulsar timing program ©. 
The position and main orbital parameters have been determined from the high-precision timing 
enabled by the strong, narrow pulses of A (Table 1). For B, an offset of 180° was added to the 
longitude of periastron and the only fitted parameters were the pulsar rotational period, P, its 
first derivative, P, and the projected semi-major axis, x B = a B sin i/c, where i is the orbital 
inclination and c is the speed of light (Table 1). 

Gravitational Physics. With their strong gravitational fields and rapid motions, DNS binaries 
exhibit large relativistic effects. General relativity and other theories of gravity can be tested 
when a number of relativistic corrections, the so-called post-Keplerian (hereafter PK) parame- 
ters, to the classical Keplerian description can be measured. In this formalism, for point masses 
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with negligible spin contributions, the PK parameters in each theory should only be functions 
of the a priori unknown neutron star masses and the classical Keplerian parameters. With the 
two masses as the only free parameters, the measurement of three or more PK parameters over- 
constrains the system, and thereby provides a test-ground for theories of gravity 0. In a theory 
that describes the binary system correctly, the PK parameters define lines in a mass-mass dia- 
gram that all intersect in a single point. Such tests have been possible to date in only two DNS 
systems, PSR B1913+16 © and PSR B1534+12 (9) (see also (123). For PSR B1913+16, the 
relativistic periastron advance, Co, the orbital decay due to gravitational wave damping, P b , and 
the gravitational redshift/time dilation parameter, 7, have been measured, providing a total of 
three PK parameters. For PSR B1534+12, Shapiro delay (|77T) . caused by passage of the pulses 
through the gravitational potential of the companion, is also visible, since the orbit is seen nearly 
edge-on. This results in two further PK parameters, r ("range") and s ("shape") of the Shapiro 
delay. However, the observed value of Pb requires correction for kinematic effects (|i2|) . so that 
PSR B 1534+ 12 provides four PK parameters usable for precise tests (9). 

Extending and improving the timing solution for A ©, we have measured A's Co and 7 and 
have also detected the Shapiro delay in the pulse arrival times of A due to the gravitational field 
of B (Fig. 2). This provides four measured PK parameters, resulting in a m A — m B plot (Fig. 3) 
through which we can test the predictions of general relativity &13W1 1 ). However, the detection 
of B as a pulsar opens up opportunities to go beyond what is possible with previously known 
DNS binary systems. Firstly, we can exclude all regions in the m^—m B plane that are forbidden 
by the individual mass functions of A and B due to the requirement smi < 1. Secondly, 
with a measurement of the projected semi-major axes of the orbits of A and B (Table 1), we 
obtain a precise measurement of the mass ratio, R{mx-,rn-Q) = mx/m-Q = xb/xa, providing 
a further constraint in the -mg plot (Fig. 3). This relation is valid for any theory of 
gravity (ITII7?t. Most importantly, the P-line is independent of strong-field (self-field) effects, 
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providing a stringent and new constraint for tests of gravitational theories il3\\ll\U5[ . 

With four PK parameters already available for tests, this additional constraint makes this 
system the most over-determined DNS binary to date and a truly unique laboratory for rela- 
tivistic gravity. Moreover, with a significant measurement of P h expected within the next few 
months, an additional PK parameter will become available. This may provide a sixth constraint 
if kinematic effects are negligible or can be isolated by proper motion and distance measure- 
ments. The position of the allowed region in Fig. 3 also determines the inclination of the orbit 
to the line-of-sight. It turns out that the system is observed nearly edge-on with an inclination 
angle i of about 87° (Table 1). 

Due to the curvature of space-time near massive objects, the spin axes of both pulsars will 
precess about the total angular momentum vector, changing the orientation of the pulsars as seen 
from Earth (176b. With the measured system parameters (Table 1), general relativity predicts 
periods of such geodetic precession of only 75 yr for A and 71 yr for B (fTTb. Hence, the relative 
orientation of the pulsars' spin axes within the system geometry is expected to change on short 
time-scales. This should lead to measurable changes in the profiles of A and B (cf. M8l\19\[ZUb ). 
and perhaps also to measurable changes in the aberration effects due to the rotation of A and B. 
Hence, we can expect that additional PK parameters can be measured that are too small to be 
detected in other binary pulsars. One example of these are the aberration terms A, B and 5 r of 
the DD timing model that we used (liilliil) . assumed to be zero for the present analysis. 

In contrast to previous tests of general relativity, we may soon have to use higher order terms 
than (f/c) 2 to describe the system accurately, the first time this has been necessary for binary 
pulsar systems. The future high precision of the measurement of Co may demand this first for 
comparison of the observed value of this parameter with theories of gravity (EZJ- Deviations 
from the value predicted by general relativity may be caused by contributions from spin-orbit 
coupling (I77I) . which is about an order of magnitude larger than for PSR B 1913+16. This 
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potentially will allow us to measure the moment of inertia of a neutron star for the first time 

Origin and evolution of the double-pulsar system. The existence of DNS binaries can be 
understood by a binary evolution scenario which starts with two main sequence stars (e.g. (IZ7V). 
The initially more massive star evolves first and eventually explodes in a supernova to form 
a neutron star. Under favorable conditions, this neutron star remains bound to its companion 
and spins down as a normal pulsar for the next 10 6 -10 7 yr. At some later time, the remaining 
(secondary) star comes to the end of its main-sequence lifetime and begins a red giant phase. 
Depending on the orbital parameters of the system, the strong gravitational field of the neutron 
star attracts matter from the red giant, forming an accretion disk and making the system visible 
as an X-ray binary. The accretion of matter transfers orbital angular momentum to the neutron 
star, spinning it up to short periods and dramatically reducing its magnetic field (|221E21l- A lim- 
iting spin period is reached due to equilibrium between the magnetic pressure of the accreting 
neutron star and the ram pressure of the infalling matter (123IL26I) . 

A crucial phase in the evolution of close DNS binaries like J0737— 3039 is the dramatic 
reduction in orbital separation that occurs when matter from the secondary star is expelled 
from the system, resulting in a very compact system consisting of a helium star and a neutron 
star (12 71 1281) . A sufficiently massive helium star will ultimately undergo a supernova explosion 
forming a young, second neutron star. If the stars remain bound following this explosion, the 
resulting system is a pair of neutron stars in an eccentric orbit with very different magnetic 
field strengths and hence spin-down properties, as in fact observed here. With masses of = 
1.34 M Q and m B = 1.25 M , A is typical of other neutron stars with measured masses i29\\10l . 
while B has a significantly smaller mass than any other. 

The time since the second supernova explosion can be estimated by comparing our mea- 
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surements of P and P for A and B, which can be used to compute their characteristic ages 
t = Pj (2P). If characteristic ages are good indicators of pulsars' true ages, we expect ta = tr, 
but the observed values are ta — 4r B . This discrepancy can be reconciled by questioning one 
or more of the assumptions inherent in the use of characteristic ages as estimates of true ages: 
a negligible birth spin period and a non-decaying magnetic dipole braking torque. At the very 
least, the post-accretion spin period of A cannot have been negligible due to details of the 
accretion process discussed above. Simple models assuming constant magnetic dipole spin- 
down (l2ol predict a post-accretion spin period for A in the range 10-18 ms, enough to explain 
the observed discrepancy in the characteristic ages. 

Probing pulsar magnetospheres. The separation of the two pulsars in their orbits (Fig. 4 
top) is typically ^900,000 km or 3 lt-s (the distance light travels in three seconds). The large 
orbital inclination means that, at conjunction, the line-of-sight to one pulsar passes within about 
0.15 lt-s of the other (Fig. 4 bottom). This is substantially smaller than the 0.45 lt-s radius of 
B's light cylinder, the point at which the co-rotation speed equals the speed of light, although 
much greater than the 0.004 lt-s light cylinder radius of A. As the pulsars move in their orbits, 
the line-of-sight from A passes through, and sweeps across, the magnetosphere of B, provid- 
ing the opportunity to probe its physical conditions. The determination of changes in the radio 
transmission properties, including the dispersion and rotation measures, will potentially allow 
the plasma density and magnetic field structure to be probed. Additionally, the ~70-yr period 
geodetic precession will cause the line-of-sight to sample different trajectories through the mag- 
netosphere of B. To a lesser extent, the 21-yr orbital precession will also cause changes in the 
trajectory. 

Close inspection of the flux density of the pulses from A (Fig. 5) reveals that a short occulta- 
tion occurs, centered upon its superior conjunction (when the line-of-sight passes only 0.15 lt-s 
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from B). The duration of the occultation is about 20-30 sec. Because the relative transverse ve- 
locity is about 660 km/s, the eclipsing region has a lateral extent of ~ 15,000 km or ~0.05 lt-s, 
about 10% of the light cylinder radius of B. As far as we can tell at present, the eclipse occurs 
in every orbit and its extent is essentially the same at 680 MHz and 1390 MHz, limiting the 
possible interpretations. One clue may come from the fact that the rate of spin-down energy 
loss from A is ~3600 times greater than that from B. In fact, at the light cylinder radius of B, 
the energy density of the relativistic wind from A is about two orders of magnitude greater than 
that from B, ensuring that the wind from A will penetrate deep into B's magneto sphere. We 
find that the energy densities due to the spin-down luminosity (assumed isotropic) emitted by 
A and the magnetic field of B are in balance at a distance ~0.2 lt-s from B, about 40% of the 
light cylinder radius of B, assuming that the magnetic field can be calculated using the standard 
dipole formula (TJDt. Within this distance, the energy field of B will dominate, while A will 
dominate outside. Even though this picture is roughly consistent with the lateral extent of the 
observed eclipse of A, it is likely that the strong departure of B's magnetosphere from an ideal 
case makes such calculations uncertain. 

Because the point of pressure-balance is deep within the magnetosphere of B, the actual 
penetration of the wind from A into B's magnetosphere will be a function of the orientation 
of the rotation and/or magnetic axes of B relative to the direction of the wind and hence will 
depend on the precessional and orbital phases of B. This is the most likely explanation for 
the large flux density and pulse- shape changes of B which are seen to vary with orbital phase 
(Fig. 1). We note that these changes are essentially the same between 680 MHz and 3030 MHz. 
It seems most likely to us that such a broadband modulation arises from the impact of energy 
in the form of particles, gamma-rays or 44-Hz electromagnetic radiation from the millisecond 
pulsar A upon the magnetosphere of B. We note that there is a strong indication of significant 
unpulsed radio emission from the system: the time-averaged pulsed flux density of the pulsars 
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is about 1.8 mJy at 1390 MHz (Table 1), compared with a total flux density at this frequency 
of 7 mJy ( J). The ~5 mJy unpulsed emission probably arises in the impact region described 
above. We find it remarkable, with much of the magnetosphere of B blown away by the wind 
of A, that B still works as a pulsar. This suggests that the radio emission is probably generated 
close to the neutron star, providing a direct constraint on the emission height. 

Conclusion. We have detected the binary companion of the millisecond pulsar J0737— 3039 
as a pulsar, making this the first known double-pulsar system. This discovery confirms the neu- 
tron star nature of the companions to recycled pulsars in eccentric binary systems and validates 
the suggested evolutionary sequences in which a companion star, having spun up the pulsar, 
forms a young pulsar in a supernova explosion 071) . The highly-relativistic nature of this com- 
pact system opens up opportunities for much more stringent tests of relativistic gravitation than 
have been possible previously. Not only have we already measured four quantities attributable 
to, and consistent with, general relativity, but the mass-ratio R is a new high-precision constraint 
that is independent of gravitational theories. Within a year or so we expect to measure the orbital 
decay due to emission of gravitational radiation. If the intrinsic value due to gravitational-wave 
damping can be extracted, it will allow tests of radiative aspects of gravitational theories mixed 
with strong-field effects. On somewhat longer time-scales of a few years we expect to detect 
several other relativistic effects such as geodetic precession of the pulsars' spin axes, spin-orbit 
coupling and other deviations, making this a superb testbed for relativity. 

The detection of the companion as a pulsar also opens up the possibility of using each pulsar 
to probe the magnetosphere of the other. The energy flux from the millisecond pulsar is strongly 
affecting the pulse emission process in the companion and eclipses of the millisecond pulsar by 
the companion are also seen. Future measurements of orbital variations in pulse shapes, ampli- 
tudes, polarization and timing over a range of radio frequencies will give fascinating insights 
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into magneto spheric processes in pulsars. 
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Table 1: Observed and derived parameters of PSRs J0737— 3039A and B using the DD timing 
model M3H111 . Standard (1 o) errors are given in parentheses after the values and are in units 
of the least significant digit(s). The parameters A, B and 8 r in the DD model were assumed to 
be zero in the analysis. The distance is estimated from the dispersion measure and a model for 
the interstellar free electron distribution (U2"l) . 



Pulsar 


PSR J0737-3039A 


PSR J0737— 3039B 


Pulse period P (ms) 


22.69937855615(6) 


2773.4607474(4) 


Period derivative P 


1.74(5) x 1(T 18 


0.88(13) x 1(T 15 


Epoch of period (MJD) 


52870.0 


52870.0 


Right ascension (J2000) 


07 h 37 m 51 s .247(2) 




Declination (J2000) 


-30°39'40".74(3) 




Dispersion measure DM (cm~ 3 pc) 


48.914(2) 


48.7(2) 


Orbital period P h (day) 


0.102251563(1) 




Eccentricity e 


0.087779(5) 




Epoch of periastron T (MJD) 


52870.0120589(6) 




Longitude of periastron uj (deg) 


73.805(3) 


73.805 + 180.0 


Projected semi-major axis x = asini/c (sec) 


1.41504(2) 


1.513(4) 


Advance of periastron Co (deg/yr) 


16.90(1) 




Gravitational redshift parameter 7 (ms) 


0.38(5) 




Shapiro delay parameter s 


0.9995(-32,+4) 




Shapiro delay parameter r (/xs) 


5.6(-12,+18) 




RMS timing residual (fis) 


27 


2660 


Flux density at 1390 MHz (mJy) 


1.6(3) 


0-1.3(3) 



Characteristic age r (My) 

Surface magnetic field strength B (Gauss) 

Spin-down luminosity E (erg/s) 

Mass function (M ) 

Distance (kpc) 

Total system mass + m B (M ) 
Mass ratio R = m^jm^ 
Orbital inclination from Shapiro s (deg) 
Orbital inclination from (R, Co) (deg) 
Stellar mass from (R, Co) (M ) 



210 

6.3 x 10 9 
5800 x 10 30 
0.29097(1) 



50 

1.6 x 10 12 
1.6 x 10 30 
0.356(3) 



-0.6 



2.588(3) 
1.069(6) 
87(3) 
87.7(-29,+17) 
1.337(5) 1.250(5) 
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Fig. 1. The intensity of radiation of B as a function of orbital longitude relative to the ascending 
node of its orbit. Each panel shows a grey -plot of intensity over a phase range of 0.1 of pulsar 
rotational period, centered on the pulsed emission. Panels (A), (B) and (C) present observations 
at 680 MHz, 1390 MHz and 3030 MHz respectively (summed over 6 orbits at 1390 MHz and 
2 orbits at each of the other frequencies), showing the similarity of the intensity variation over 
the wide frequency range and the changing pulse-shape with orbital longitude. Longitude in 
this diagram (and in Figs. 2 and 5) is strictly the sum of the longitude of periastron and the true 
anomaly. The vertical line at longitude 270° represents the longitude of inferior conjunction of 
B, when the two stars are at their closest on the sky and B is the nearer to the Earth. 
Fig. 2. The effect of the Shapiro delay caused by the gravitational potential of B seen in the 
timing residuals of A. Residuals were averaged into 75 equal bins of orbital longitude. Top: 
timing residuals obtained by subtracting the model defined in Table 1 from the observed TOAs. 
Bottom: as top display, but with the Shapiro delay parameters r and s set to zero. 
Fig. 3. The observational constraints upon the masses m A and m B . The colored regions are 
those which are excluded by the Keplerian mass functions of the two pulsars. Further con- 
straints are shown as pairs of lines enclosing permitted regions as predicted by general relativity 
(see text for details): (a) the measurement of the advance of periastron to, giving the total mass 
mA+m B = 2.588±0.003 M Q (dashed line); (b) the measurement of R = m A /m B = xb/xa = 
1.069 ±0.006 (solid line); (c) the measurement of the gravitational redshift/time dilation param- 
eter 7 (dot-dash line); (d) the measurement of Shapiro parameter r giving m B = 1.2 ± 0.3 M 
(dot-dot-dot-dash line) and (e) Shapiro parameter s (dotted line). Inset is an enlarged view of 
the small square which encompasses the intersection of the three tightest constraints, with the 
scales increased by a factor of 16. The permitted regions are those between the pairs of parallel 
lines and we see that an area exists which is compatible with all constraints, delineated by the 
solid blue region. 
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Fig. 4. The physical configuration of the binary system, at conjunction, as on 19 August 2003 
(MJD 52870), showing the relative sizes of the two orbits and B's magneto sphere. At conjunc- 
tion, the two neutron stars are separated by ~2.8 lt-s or about 800,000 km. Top: view from 
above the orbital plane with the Earth to the right. The shaded segments indicate orbital phases 
where B is detected strongly (see Fig. 1). Also shown is the apsidal line, which is the major 
axis of the orbits, and the line-of-nodes, which is the intersection of the orbital plane and the 
plane through the center of mass of the system which is normal to the line-of-sight. Bottom: 
view from the side, showing the passage of the line-of-sight from A to the Earth through the 
magnetosphere of B. The approximate position of the pressure balance between the relativistic 
wind from A and the magnetic field of B is indicated. 

Fig. 5. The variation in flux density of A (in arbitrary units) at 680 MHz and 1390 MHz, around 
superior conjunction (i.e. longitude 90°). The data are presented with 5-sec time resolution and 
show the eclipse of the pulsar by the magnetosphere of B. 
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